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The symmetry energy for nuclear matter and its relation to the neutron skin in finite nuclei is 
discussed. The symmetry energy as a function of density obtained in a self-consistent Green function 
approach is presented and compared to the results of other recent theoretical approaches. A partial 
explanation of the linear relation between the symmetry energy and the neutron skin is proposed. 
The potential of several experimental methods to extract the neutron skin is examined. 
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(S| . I- INTRODUCTION 

The nuclear symmetry energy plays a central role in a variety of nuclear phenomena. It determines to a large extent 
Q ' the equation of state (EoS) and the proton fraction of neutron stars Q], the neutron skin in heavy nuclei 0, it enters as 
an input in the analysis of heavy ion reactions etc. Its value at nuclear saturation density, S(po — 0.17fm~ 3 ) sa 

30 MeV, seems reasonably well established, both empirically as well as theoretically; still different parametrizations 
of relativistic mean field (RMF) models (which fit observables for isospin symmetric nuclei well) lead to a relatively 
y—i • wide range of predictions for the symmetry energy, 24-40 MeV. However, predictions for its density dependence show 
^ ' a substantially larger variation. 

Recently it has been pointed out by several authors 0, 0, El that there exists a strong correlation between the 
neutron skin, AR — R n — R p , and the density derivative of the EoS of neutron matter near saturation density. 
Subsequently in a more detailed analysis in the framework of a mean field approach Furnstahl Q demonstrated that 
in heavy nuclei there exists an almost linear empirical correlation between theoretical predictions in terms of various 
(f-) \ mean field approaches to S(p) (i.e., a bulk property) and the neutron skin, AR (a property of finite nuclei). 

This observation has contributed to a renewed interest in an accurate determination of the neutron skin in neutron 
rich nuclei for several reasons. First precise experimental information on the neutron skin could help to further 
constrain interaction parameters that play a role in the calculation of the symmetry energy [tJ • Furthermore a precise 
value of the neutron skin is required as an input in several processes of physical interest, e.g. the analysis of energy 
shifts in deeply bound pionic atoms 0, and in the analysis of atomic parity violation experiments (weak charge) 0. It 



has been shown that the calculated symmetry energy is quite insensitive to details of modern realistic NN interactions 
On the other hand the symmetry energy and in particular its density dependence can vary substantially with 
the many-body approximations employed. For instance the results of lowest order BHF and variational calculations 
do not seem to agree well. 

The aim of this paper is threefold: First we address the sensitivity of the symmetry energy to various many-body 
approximations. To this end we present results of a calculation of S(p) using the self-consistent Green function (SCGF) 
approach and compare the result with various other theoretical approaches. Secondly we will try to provide some 
new insight in the origin of the "Furnstahl" relation; for this purpose we use the Landau-Migdal effective interaction 
in the mean field approximation. Finally in view of the large variety of existing and proposed experimental methods 
to determine the neutron skin AR we examine the merits of some recently proposed methods that seem to be of 
potential interest to provide more accurate information on the neutron skin in the near future. 

Section ^ is devoted to an overview of theoretical approaches to the symmetry energy, and a new calculation in 
terms of the self-consistent Green function approach is presented. In Section IlIII an interpretation of the Furnstahl 
relation is presented in terms of the Landau-Migdal approach. Section llVl contains an overview of various experimental 
methods to deduce information on the neutron skin and Section^contains a short discussion of implications for other 
physical processes where the information on the neutron skin is required as an input. 



II. THE SYMMETRY ENERGY IN NUCLEAR MATTER 



The symmetry energy S(p) is defined in terms of a Taylor series expansion of the energy per particle for nuclear 
matter in terms of the asymmetry a — (N — Z)/A (or equivalently the proton fraction x — Z/A), 

E(p, a) = E(p, 0) + S(p)a 2 + 0(a 4 ) + . . . . (1) 
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It has been shown jl2lll3| that deviations from the parabolic law in Eq.QJ, i.e. terms in a 4 , are quite small. 

Near the saturation density p the energy of isospin-symmetric matter E(p,0) and the symmetry energy can be 
expanded as 

E{p,0) = £o + -^7 (p-po) 2 + -.. , (2) 
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and 
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S (p) = o — ETT^ — l<*=o = fl 4 + —(P-Po) + t^-jC/'-Po) 2 + ■ ■ •• (3) 



The parameter a4 is the symmetry energy at equilibrium and the slope parameter p governs the density dependence. 

The relevance of the nuclear matter results in part depends on the question whether there is a surface contribution 
to the symmetry energy for finite nuclei. In ref.|l4j it was found that the latter is of minor importance, which has 
also been confirmed in [2|. 

A. Self-consistent Green function and Brueckner approach 

In this section we describe the calculation of the symmetry energy in the SCGF approach. Since the latter can 
be considered as a generalization of the lowest order Brueckner-Hartree-Fock (BHF) method we start with a brief 
discussion of the symmetry energy in the latter case. 

1. Symmetry energy in BHF 

In the Brueckner-Hartree-Fock approximation, the Brueckner-Bethe-Goldstone (BBG) hole-line expansion is trun- 
cated at the two hole-line level. The short-range NN repulsion is treated by a resummation of the particle-particle 
ladder diagrams into an effective interaction or G-matrix. Self-consistency is required at the level of the BHF single- 
particle spectrum e BHF (k), 

e BHF {k) = jf_ + y Re< kk , lG{uJ = e BHF {k) + £ BHF mkk , > (4) 

2m 

k'<k F 

In the standard choice BHF the self-consistency requirement (0} is restricted to hole states (k < k F ) only, while the 
free spectrum is kept for particle states k > kp. The resulting gap in the s.p. spectrum at k = kp is avoided in 
the continuous-choice BHF (ccBHF), where Eq. J1J is used for both hole and particle states. The continuous choice 
for the s.p. spectrum is closer in spirit to many-body Green's function perturbation theory. Moreover, recent results 
indicate |15l Ha ] that the contribution of higher-order terms in the hole-line expansion is considerably smaller if the 
continuous choice is used. 

The BHF energy per nucleon can be easily evaluated for both symmetric nuclear matter (SNM) and pure neutron 
matter (PNM) using the energy sum-rule, 

E '' ' ''''''' 1 +e BHF (k))e(k F -k), (5) 



A pj (2tt) 3 \2m 

where p = ^f- is the density and d is the isospin degeneracy (d = 1(2) for PNM (SNM)). The symmetry energy S(p) 
is obtained as the difference between PNM and SNM energies for the same density. 

We first performed ccBHF calculations with the Reid93 interaction, including partial waves with J < 4 in the 
calculation of the G-matrix. The results are presented in Fig^ where the symmetry energy S is decomposed into 
various contributions as suggested in ref.|l2| and shown as a function of nucleon density p. The contribution 5km of 
the kinetic energy to the BHF symmetry energy is given by the free Fermi-gas expression |44| 

3 22/ 1 \ 

Skin = Ekin,PNM ~ Ekin,SNM = —Z (37T ) 3 p 3 1 ) , (6) 

10m \ 23/ 

and it determines to a large extent the density dependence of S. In Fig^ we also show the symmetry potential 
Spot — S — Skin, which is much flatter, and the contributions to S po t from both the isoscalar (T = 0) and isovector 
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FIG. 1: The symmetry energy S (full line) and the contributions to S from the kinetic (dash-dotted line) and potential energy 
(dashed line), calculated within a ccBHF scheme and using the Reid93 interaction. Also shown (dotted lines) are the T — 
and T = 1 components of the potential energy contribution. 

(T = 1) components of the interaction. Over the considered density range S po t is dominated by the positive T = 
part. The T = partial waves, containing the tensor force in the 3 S±- 3 Di channel which gives a major contribution 
to the potential energy in SNM, do not contribute to the PNM energy. The T = contribution peaks at p m 0.3 
fm — . The decrease of this contribution at higher densities is compensated by the increase of the T — 1 potential 
energy, with as a net result a much weaker density-dependence of the total potential energy. (It should be noted that 
inclusion of a three- nucleon interaction in general leads to a substantial increase for the slope parameter po [l3j|.) 

2. Symmetry energy in SCGF approach 

In recent years several groups have considered the replacement of the BBG hole-line expansion with self-consistent 
Green's function (SCGF) theory [H |U El |H |H . In [2^, 0| the binding energy for SNM was calculated within 
the SCGF framework and using the Reid93 potential. In the present paper we have extended these calculations to 
PNM and calculated the corresponding symmetry energy. Details of a technical nature can be found in [2^| . 

A SCGF calculation differs in two important ways from a BHF- calculation. Firstly, within SCGF particles and 
holes are treated on an equal footing, whereas in BHF only intermediate particle (fc > /cp) states are included in 
the ladder diagrams. This aspect ensures thermodynamic consistency, e.g. the Fermi energy or chemical potential of 
the nucleons equals the binding energy at saturation (i.e., it fulfills the Hugenholz-van Hove theorem). In the low- 
density limit BHF and SCGF coincide. As the density increases the phase space for hole-hole propagation is no longer 
negligible, and this leads to an important repulsive effect on the total energy. Secondly, the SCGF generates realistic 
spectral functions, which are used to evaluate the effective interaction and the corresponding nucleon self-energy. The 
spectral functions include a depletion of the quasi-particle peak and the appearance of single-particle strength at 
large values of energy and momentum, in agreement with experimental information from (e, e'p) reactions. This is in 
contrast with the BHF approach where all s.p. strength is concentrated at the BHF-energy as determined from Eq. 

In a SCGF approach the particle states (fc > hp), which are absent in the BHF energy sum rule of Eq. (J5J, do 
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FIG. 2: The total-energy per particle for symmetric nuclear matter (left panel) and pure neutron matter (central panel) for 
the Reid93-interaction. The dashed line refers to a ccBHF calculation, the full line to a SCGF calculation. The right panel 
displays the symmetry energy in these two approaches. 



contribute according to the energy sum rule 

E d f d 3 k [ SF , / k 2 \ . 

expressed in terms of the nucleon spectral function Sh(k,u>). 

The results for the ccBHF and SCGF calculations for both SNM and PNM are compared in the left and central 
panels of Fig[21for the Reid93 interaction. The inclusion of high-momentum nucleons leads roughly to a doubling of 
the kinetic and potential energy in SNM, as compared to BHF. As seen in Fig|3 the net result for the total energy 
of SNM is a repulsive effect, increasing with density (2j|. This leads to a stiffer equation of state, and a shift of the 
SNM saturation density towards lower densities. The above effects are dominated by the tensor force (the isoscalar 
3 Si — 3 D\ partial wave). Consequently, the effects are much smaller in PNM. 

The corresponding symmetry energy, shown in the right panel of Fig|2 is dominated by the shift in the total energy 
for SNM, and lies below the ccBHF symmetry energy in the entire density-range. At po — 0.16fm~ 3 the symmetry 
energy parameter C14 is reduced from 28.9 MeV to 24.9 MeV, while the slope po remains almost the same (from 2.11 
MeVfm" 3 to 1.99 MeVim" 3 ). 



B. Comparison of symmetry energy in other approaches 

1. Calculations with realistic NN forces 

Engvik et al. have performed lowest-order Brueckner-Hartree-Fock calculations in SNM and PNM for all 
"modern" potentials (CD-Bonn, Argonne vl8, Reid93, Nijmegen I and II), which fit the Nijmegen NN scattering 
database with high accuracy. They concluded that for small and normal densities the symmetry energy is largely 
independent of the interaction used, e.g. at po the values of S vary around an average value of 04=29.83 MeV by 
about 1 MeV. At larger densities the spread becomes larger; however, the symmetry energy keeps increasing with 
density, in contrast to some of the older potentials like Argonne vl4 and the original Reid interaction (Reid68) for 
which S(p) tended to saturate at densities larger than p = 0.4 fm~ 3 . Some insight into the microscopic origin of the 
symmetry potential was obtained by Zuo et al. ^3] w ho decomposed the symmetry energy into contributions from 
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TABLE I: Results for the symmetry energy parameters 04 and po from the variational calculations of Refs.|l7 | .|l8 j| using the 
Argonne and Urbana NN potentials, in combination with Urbana models for the three-nucleon interaction. The last column 
includes a relativistic boost correction Sv and the adjusted UIX* three-nucleon interaction. 





Avl4 


Avl4+UVIII 


Uvl4 


Uvl4+UVIII 


Avl8 


Avl8+UIX 


Av18+<5?j+UIX* 


a 4 [MeV] 


24.90 


27.49 


26.39 


28.76 


26.92 


29.23 


30.1 


po [MeV fm~ 3 ] 


2.02 


2.71 


2.38 


3.04 


1.95 


3.24 


2.95 



kinetic and potential energy. The BHF calculations in 12] used the Argonne vl4 and the separable Paris interaction. 
In FigQwe showed that the use of the modern Reid93 potential leads to essentially the same conclusions. 

Detailed studies for SNM and PNM using variational chain summation techniques were performed by Wiringa et 
al. ^3 f° r the Argonne Avl4 and Urbana Uvl4 NN interaction, in combination with the Urbana UVIII three-nucleon 
interaction (TNI), and by Akmal et al. |l8l for the modern Avl8 NN potential in combination with the ULX-TNI. 
Results for 04 and po, extracted from pjj and [lj|, are shown in Tabled The inclusion of TNI stiffens the EoS for 
both SNM and PNM, and increases considerably the value of the symmetry energy 0,4 and its slope po at the empirical 
density po- The effect of including a relativistic boost correction 5v (in combination with a refitted TNI) on the values 
of 04 and po is sizeable as well. 

The symmetry energy has also been computed in the Dirac-Brueckner-Hartree-Fock (DBHF) approach l2fij. 
In relativistic approaches the symmetry energy generally is found to increase almost linearly with density, and more 
rapidly than in the non-relativistic case. This difference can be attributed to two effects. First the covariant kinetic 
energy which is inversely proportional to •v/fc^ + to* 2 is larger because of the decreasing Dirac mass, to*, with in- 
creasing density. Secondly the contribution from rho-exchange appears to be larger than in the non-relativistic case 




2. Mean-field approach using effective interactions 



Since the Furnstahl relation has been verified mainly in terms of mean-field models we discuss some results obtained 
in these approaches, which in general are based upon a parameterized effective interaction. 

Brown [5( has investigated proton and neutron radii in terms of the non-relativistic Skyrme Hartree-Fock (SHF) 
model. First he noted that a certain combination of parameters in the SHF is not determined well by a fit to ground 
state binding energies, and that a wide range of predictions for the EoS for PNM is obtained. He also pointed out 
a direct correlation between the derivative of the neutron matter EoS (i.e., basically the symmetry energy coefficient 
po) and the neutron skin in 208 Pb. 

Covariant approaches are in general based upon either a covariant lagrangian with a, lo and p exchange (and 
possibly other mesons) 0, , or on the use of contact interactions , solved as an energy density functional (EDF) 
in the Hartree-Fock approximation. Sets of model parameters are determined by fitting bound state properties of 
nuclei. Specifically the isovector degree of freedom is determined by the exchange of isovector mesons; in case of 
rho-meson exchange the (positive definite) contribution to S is given by 

k 2 a 2 

J? 2 rr| 

and its potential energy contribution to po, which scales with that for 04, is g^Aj- [2j. Typical values obtained for po are 

around 4-6 MeV fm -3 , and 04 ~ 30-36 MeV, i.e. considerably larger than in non-relativistic approaches (a large part 
of the enhancement can be ascribed to the fact that the kinetic contribution is larger, because to* < to). Recently in 
[28l |2^ | this approach was extended by inclusion of the isovector-scalar partner, S, of the isoscalar scalar a— meson. 

Because of the presence of the Lorentz factor m* /E in the scalar potential contribution, ~ — gs^gf , which decreases 

with increasing density its inclusion leads to an even larger net value for po psl l29j . 

Since explicit pion exchange is usually not included in the mean field approaches it is difficult to make a meaningful 
comparison with microscopic ones. In fact it can be argued that in contrast to isoscalar properties the long-range 
pion exchange could play an essential role in determining the isovector properties 0- 
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FIG. 3: Overview of several theoretical predictions for the symmetry energy S: Brueckner-Hartree-Fock (continuous choice) 
with Reid93 potential (circles), self-consistent Green function theory with Reid93 potential (full line), variational calculation 
from |r?| with Argonne Avl4 potential (dashed line), Dirac-Brueckner-Hartree-Fock calculation from |25l| (triangles), relativistic 
mean-field model from |28j| (squares), effective field theory from |3(J (dash-dotted line). 



3. Effective field theory 

Recently the density dependence of the symmetry energy has been computed in chiral perturbation effective field 
theory, described by pions plus one cutoff parameter, A, to simulate the short distance behavior 30] . The nuclear 
matter calculations have been performed up to three-loop order; the resulting EoS is expressed as an expansion in 
powers of kp, and the value of A » 0.65 GeV is adjusted to the empirical binding energy per nucleon. The value 
obtained in this approach for 04 = 33 MeV is in remarkable agreement with the empirical one; at higher densities 
(p > 0.2 fm~ 3 ) a downward bending is predicted. However, in its present form the validity of this approach is clearly 
confined to relatively small values of the Fermi momentum, i.e. rather low densities. It is interesting to note that there 
are relatively small (large) contributions to 04 coming from one-pion exchange Fock diagram (three-loop diagrams 
with either two or three medium insertions). 

4- Comparison 

To summarize the present status in Fig. ljll B 4f) various results of the approaches for S(p) discussed above are 
compared. As noted above one sees that the covariant models yield a much larger increase of S with the density 
than the non-relativistic approaches. The LOBHF leads to a higher value of S than both variational and the SCGF 
method which include more correlations; that the SCGF result is close to the variational approach may be fortuitous. 
Effects from three-body forces are not included. 
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III. RELATIONSHIP BETWEEN SYMMETRY ENERGY AND AR 



Brown [5| and Furnstahl |2J have pointed out that within the framework of mean field models there exists an 
almost linear empirical correlation between theoretical predictions for both a 4 and its density dependence, po? and 
the neutron skin, AR = R n — R p , in heavy nuclei. This observation suggests an intriguing relationship between a 
bulk property of infinite nuclear matter and a surface property of finite systems. 

Here wc wish to address this question from a different point of view, namely in the spirit of Landau-Migdal approach. 
Let us consider a simple mean- field model (see, e.g., |31| ) with the Hamiltonian consisting of the single-particle mean 
field part H and the residual particle- hole interaction H p -h- 

H = H + H P -h, H = J2(T a + U(x a )), (9) 

a 

U(x) = U (x) + U 1 (x) + U c (x), (10) 
U (x) = U {r) + U so (x); Ui(x) = ±S pot (r)r( 3 ) ; U c {x) = \U c (r){\ - t^), x = (r,a,r) 

Here, the mean field potential U(x) includes the phenomenological isoscalar part Uq(x) along with the isovector U\(x) 
and the Coulomb Uc{x) parts calculated consistently in the Hartree approximation; Uo(r) and U so {x) = U so (r)a ■ I 
are the central and spin-orbit parts of the isoscalar mean field, respectively; S po t(r) is the symmetry potential (the 
potential part of the symmetry energy). 

In the Landau-Migdal approach the effective isovector particle-hole interaction H p -h is given by 

Hph = J2( F ' + G'a a a b )f a f b 5(r a - f b ), (11) 

a>b 

where F' and G' are the phenomenological Landau-Migdal parameters. 

The model Hamiltonian H in Ea. (|l(Jfl preserves isospin symmetry if the condition 

[#,T<->] = £><->, (12) 

is fulfilled, where = £ , U [ J ] = ^ U c (r a )Tt ] . With the use of Eqs. JTUJl , dJ , JUJl the condition ea.lfT^I in 

a a 

the RPA formalism leads to a selfconsistency relation between the symmetry potential and the Landau parameter F' 

M 

S pot {r)=2F'n(-\r), (13) 

where n(~>(r) = n n (r) — n p (r) is the neutron excess density. Thus, in this model the depth of the symmetry potential 
is controlled by the Landau-Migdal parameter F' ( analogously to the role played by the parameter g 2 p in relativistic 
mean field models). 

Spotif) i s obtained from Ea. (|13|l by an iterative procedure; the resulting dependence of AR on the dimcnsionlcss 
parameter /' = F'/(300 MeVfm 3 ) shown in fig. 0] indeed illustrates that AR depends almost linearly on /'. Then 
with the use of the Migdal relation |33J which relates symmetry energy and /', 

04 = ^(1 + 2/'), (14) 

a similar almost linear correlation between the symmetry energy, 0,4,, and the neutron skin is obtained. 

To get more insight in the role of /' we consider small variations SF'. Neglecting the variation of n^(r) with 
respect to SF' one has a linear variation of the symmetry potential: SS pot (r) — 2SF'n^ (r). Then in first order 
perturbation theory, such a variation of S poi causes the following variation of the ground-state wave function 

s 

with s labeling the eigenstates of the nuclear Hamiltonian and a single-particle operator — J2 n ^~H r a) T a 3 ' > ■ 

a 

Consequently the variation of the expectation value of the single-particle operator V"(~) = r a T a 3 ^ with (0|V"(~) |0) = 



NR^ — ZR p can be written as 
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FIG. 4: Neutron skin in 208 Pb vs. the Landau-Migdal parameter /'. 

In practice the sum in Eq. (|16f) is exhausted mainly by the isovector monopole resonance of which the high excita- 
tion energy (about 24 MeV in 208 Pb) justifies the perturbative consideration. We checked that Eq. is able to 
reproduce directly calculated S(AR) shown in fig.0|with the accuracy of about 10%. As a result a simple microscopic 
interpretation of the linear correlation between the neutron skin thickness and Landau parameter F' is obtained. 



IV. EXPERIMENTAL METHODS TO DETERMINE AR 



A variety of experimental approaches have been explored in the past to obtain information on AR. To a certain 
extent all analysis contain a certain model dependence, which is difficult to estimate quantitatively. It is not our 
intention to present a full overview of existing methods for the special case of 208 Pb. In particular the results obtained 
in the past from the analysis of elastic scattering of protons and neutrons have varied depending upon specifics of 
the analysis employed. At present the most accurate value comes from a recent detailed analysis of the elastic proton 
scattering reaction at E = 0.5 — 1 GeV [4lj. and of neutron and proton scattering at E = 40 — 100 MeV 42]. For 
details we refer to these papers. Here we restrict ourselves to a discussion of the some less well known methods that 
have the potential to provide more accurate information on the neutron skin in the future. In particular the excitation 
of isovector giant resonances through the restoring force contain information about the symmetry energy. 



A. Isovector giant resonances 



We begin with a brief overview of the study of excitation of isovector giant resonances. Sum rules for the latter 
contain direct information on AR. 

1. Giant Dipole Resonance (GDR) 

In the past the excitation of the isovector giant dipole resonance (GDR) with isoscalar probes has been used to 
extract AR/R 34] . In the DWBA optical model analysis of the cross section the neutron and proton transition 
densities are needed as an input. In the Goldhaber- Teller picture 

2Ni dpi 

9^) = -n—— (17) 

with k the oscillation amplitude and (i — p,n) . Assuming ground state neutron and proton distributions of the form 
(x = (N- Z)/A) 

Pi (r) = i(l ± x T ix)p{r - c(l ± jx/3)). (18) 



9 



While for N = Z the transition density vanishes, for N > Z the isovector transition density is finite 

A t \ N_Z ( dp cd 2 p 

where 7 is related to AR, 7 = 2 ^ A Z ^ AR/R . 

Excitation of the GDR by alpha particle scattering (isoscalar probe) the corresponding transition optical potential is 
given by 

AU tr = n 1 — — (— + -— ). (19) 

By comparing the experimental cross section with the theoretical one (calculated as a function of the ratio AR/ R) 
the value of AR/R can be deduced. 

It is difficult to make a quantitative estimate of the uncertainty in the result coming from the model dependence of 
the approach. In the analysis several assumptions must be made, such as the radial shape of the density oscillations 
and about the actual values of the optical model parameters. 



2. Spin-dipole Giant Resonance 



Recently it has been proposed to utilize the excitation of the spin-dipole resonance, excited in charge exchange 
reactions, to determine the neutron skin; in fact the method has been applied to obtain information on the variation 
of the neutron skin in the Sn isotopes with isotope number |35j . For the relevant operator, ^criT^riYifa), the 
summed AL = 1 strength is 

S<-> - = C{NR 2 n - ZR 2 p ). (20) 

Here and £(+) are the spin-dipole total strengths in /?( ) and channels, respectively; C is the factor 

depending on the definition of the spin-dipole operator (in the definition of Ref. [2] C = 1/4-7T, we use here C = 1). 
Because could not be measured experimentally, the model-dependent energy-weighted sum rule was invoked in 
the analysis to eliminate S^ + \ Let put = (that seems to be a very good approximation for 208 Pb) and ask the 
question what experimental accuracy for S^' 1 is needed to determine the neutron skin to a given accuracy. With 

= (N- Z)R 2 + 2NR P AR (21) 
the ratio of the second term on the rhs to the first one in case of 208 Pb is 



2NAR/({N ~ Z)Rp) « 5.7AR/R P . 

Therefore, for R p = 5.5 fm and AR = 0.2 fm the second term is only 25% of the first one and one needs 3% accuracy 
in £( - ) to determine AR with 10% accuracy. Because the SDR strength is spread out and probably has a considerable 
strength at low-energy the results for the AR can be only considered as qualitative with a relatively large uncertainty 
(of order 30-50%). 



3. Isobaric analogue state 

The dominant contribution to the energy weighted sum rule (EWSR) for Fermi excitations by the operator T(~> 
comes from the Coulomb mean field [3q 

(EWSR) F = J U c (r)n^(r)d s r, (22) 

The Coulomb mean field Uc{r) resembles very much that of the uniformly charged sphere, being inside a nucleus a 

Ze 2 

quadratic function: Uc{t) = (3 — (r/R c ) 2 ), r < R c . It turns out that if one extends such a quadratic dependence 

2R C 

also to the outer region r > R c (instead of proportionality to R c /r), it gives numerically just very small deviation in 
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(EWSR)f (less than 0.5%, due to the fact, that the difference and its first derivative go to at r = R c and >(r) 
is exponentially decreasing at r > R c ). Using such an approximation, one gets: 

{ EWSIt),*{N-QAc(l-^£^) (23) 

3Ze 2 

with A c = , and given in Ea. (|2"TJ) . 

2R C 

Since the IAS exhausts almost 100% of the NEWSR and EWSR, one may hope to extract from the IAS energy. 
However, the term depending on S^' contributes only about 20% to {EWSR)f, and as a result, the part of S(-> 
depending on AR contributes only about 4% to (EWSR)f (in 208 Pb). ^From the experimental side, the IAS energy 
can be determined with unprecendently high accuracy, better than 0.1%. Also, from the experimentally known charge 
density distribution the Coulomb mean field Uc{r) can be calculated rather accurately, and hence one can determine 
the small difference between Eqs.(|23l) and (|22[1 . But at the level of 1% accuracy several theoretical effects discarded in 
Ea. (|22|l come into play (see, e.g., |3g) which makes the reliability of such a method questionable. On the other hand 
in a forth coming paper we will show that for an isotopic chain the excitation of the IAS can be used as a quantitative 
tool to obtain the variation of AR with neutron number. 



B. Anti-protonic atoms 

Recently neutron density distributions were deduced from anti-protonic atoms [37] . The basic method determines 
the ratio of neutron and proton distributions at large differences by means of a measurement of the annihilation prod- 
ucts which indicates whether the antiproton was captured on a neutron or a proton. In the analysis two assumptions 
are made. First a best fit value for the ratio Ri of the imaginary parts of the free space pp and pn scattering lengths 
equal to unity is adopted. Second in order to reduce the density ratio at the annihilation side to a a ratio of rms radii 
a two-parameter Fermi distribution is assumed. The model dependence introduced by this assumptions is difficult to 
judge. Since a large number of nuclei have been measured one may argue that the value of Rj is fixed empirically. 

C. Parity violating electron scattering 

Recently it has been proposed to use the (parity violating) weak interaction to probe the neutron distribution. This 
is probably the least model dependent approach |3S| . The weak potential between electron and a nucleus 

V(r) = V(r)+ l5 A(r), (24) 

where the axial potential A(r) = ■^f^pwi 7 ')- The weak charge is mainly determined by neutrons 

pw(r) = (1 - 4sin 2 6 w )p P (r) - Pn(r), (25) 

with sin 2 6w « 0.23. In a scattering experiment using polarised electrons one can determine the cross section asym- 
metry |38| which comes from the interference between the A and V contributions. Using the measured neutron form 
factor at small finite value of Q 2 and the existing information on the charge distribution one can uniquely extract 
the neutron skin. Some slight model dependence comes from the need to assume a certain radial dependence for the 
neutron density, to extract R n from a finite Q 2 form factor. 

V. DISCUSSION OF AR FOR 208 PB AND SOME IMPLICATIONS 

In table [V] we present a summary of some recent results on AR in 208 Pb. One sees that (with the exception of the 
analysis of proton and neutron pickup reactions in terms of mean field orbitals in |39| | ) all recent results are consistent 
with AR ~ 0.13 ± 0.03 fm. Therefore it appears that the data agree with the result of conventional Skyrme model 
approach but seem to disagree with the results of the RMF models considered in . On the basis of the correlation 
plot between AR and po shown in |2( on would then conclude that a small value for po ~ 2MeV/fm 3 is preferred over 
larger values prdedicted in RNF approaches. 

In several processes of physical interest knowledge of AR plays a crucial role and in fact a more accurate value 
could lead to more stringent tests: 
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TABLE II: Summary of recent results for AR in 20S Pb 



method 


AR [fm] error [fm 


ref 


giant dipole resonance excitation 


0.19 


0.09 


[34, 40] 


neutron/proton pickup 


0.51 




f39] 


{p,p') at 0.5-1.04 GeV 


0.097 


0.014 


[41] 


nucleon scattering (40-200 MeV) 


0.17 




[42] 


anti-protonic atoms 


0.15 


0.02 


f37] 


parity violating electron scattering 


planned 


0.05 


[38] 



(i) The pion polarization operator (the s-wave optical potential) in a heavy nucleus H(uj, p p ,p n ) — ~T + (uj)p — 
T~(Lu)(p n — pp) has mainly an isovector character (T + (m 7r ) ~ 0). Parameterizing the densities by Fermi shapes for 
the case of 208 Pb the main nuclear model dependence in the analysis comes from the uncertainty in the value of AR 
multiplying T~ . 

(ii) Parity violation in atoms is dominated by Z— boson exchange between the electrons and the neutrons |lfl l43| . 
Taking the proton distribution as a reference there is a small socalled neutron skin (ns) correction to the parity 
non-conserving amplitude, <5i?p UC , for, say, a 6si/2 — > 7s\/2 transition, which is related to AR as (independent of 
the electronic structure) 

^ (26) 

-&pnc < Hp 

In 133 Cs it amounts to a 5E/E k, —(0.1 — 0.4)% depending on whether the non-relativistic or relativistic estimates 
for AR are used The corresponding uncertainty in the weak charge Qw is —(0.2 — 0.8)cr. 

(iii) The pressure in neutron star matter can be expressed as in terms of the symmetry energy and its density 
dependence Q 

P(p,x) = = p 2 [E > {p! 1/2) + s / ( )(1 _ 2x f + ...]. (27) 

dp 

By using the beta equilibrium condition in a neutron star, p e = p n — p p = — ^ E q P ^ , and the result for the electron 
chemical potential, p e — 3/4hcx(37r 2 px) 1 ^ 3 , one finds the proton fraction near saturation density, po, to be quite 
small, xq ~ 0.04. Hence the pressure at saturation density can be approximated as 

P(p ) = Po (l - 2x )(p S'(p )(l - 2x Q ) + S( Po )x ) ~ p 2 S'(p ). (28) 

At higher densities the proton fraction increases; this increase is more rapid in case of larger po @- While for the 
pressure at higher densities contributions from other nuclear quantities like compressibility will play a role in 
it was argued that that there is a correlation of the neutron star radius and the pressure which does not depend 
on the EoS at the highest densities. Numerically the correlation can be expressed in the form of a power law, 
R M ~ C(p, ^)( Mc vfa-^ )°' 25 km ' wri ere C(p = 1.5p Q ,M = 1.4M soiar ) ~ 7. This shows that a determination of a 
neutron star radius would provide some constraint on the symmetry properties of nuclear matter. 



VI. CONCLUSION 



In this paper we have discussed the bulk symmetry energy, and compared various approaches to compute it as a 
function of density. Because the tensor interaction plays an important role the symmetry energy is sensitive to details 
of the treatment of the many-body correlations. It was shown that the self-consistent Green function approach in 
which more correlations are included than in lowest order BHF leads to a smaller value of the symmetry energy. The 
relatively large values for po obtained in the relativistic mean field approach can be associated with an effective mass 
effect. 

We showed that the phenomenological almost linear relationship between symmetry energy and neutron skin in 
finite nuclei observed in mean field calculations could be understood in terms the Landau-Migdal approach. Finally 
we compared several experimental tools of potential interest for the determination of the neutron skin. 
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